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ABSTRACT 

Power spectra of Large Magellanic Cloud (LMC) emission at 24, 70 and 160 
/xm observed with the Spitzer Space Telescope have a two-component power-law 
structure with a shallow slope of —1.6 at low wavenumber, /c, and a steep slope 
of —2.9 at high k. The break occurs at k~^ ~ 100 — 200 pc, which is interpreted 
as the line-of-sight thickness of the LMC disk. The slopes are slightly steeper for 
longer wavelengths, suggesting the cooler dust emission is smoother than the hot 
emission. The power spectrum covers ~ 3.5 orders of magnitude and the break 
in the slope is in the middle of this range on a logarithmic scale. Large-scale 
driving from galactic and extragalactic processes, including disk self-gravity, spi- 
ral waves and bars, presumably cause the low-fc structure in what is effectively a 
two-dimensional geometry. Small-scale driving from stellar processes and shocks 
cause the high-A; structure in a 3D geometry. This transition in dimensionality 
corresponds to the observed change in power spectrum slope. A companion paper 
models the observed power-law with a self-gravitating hydrodynamics simulation 
of a galaxy like the LMC. 



Subject headings: ISM: structure — galaxies: ISM — Magellanic Clouds — In- 
frared: ISM 
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Introduction 



The Large Magellanic Cloud (LMC) is the first spiral galaxy in which an elongated 
feature of stars was identified. In a remarkable drawing of the LMC as seen with the naked eye 
made by Sir John Herschel, a bar is clearly shown, which Sir John termed an "axis of light" . 
Also faithfully represented in that drawing is a prominent spiral arm to the north, as well as 
two "embryonic" arms. The LMC belongs to the de Vaucouleurs classification bin SB(s)m, 
and is ap proximately 50 kpc di stant (based on the distance modulus m — M = 18.50 ± 0.10, 
following iFreedman et al.ll200ll ). Studies of the LMC concur that the inclination to the line- 



of-sight of the (warped) disk lies in the range 30°-50° while the posit ion angle of the line 



of nod es lies between 120° and 150°; an exce l lent r eview is provided by van der Marel et al. 



((20081). The study by 



van 



der Marel fc Cionil J200lh using DENIS and 2MASS surveys yields 
an inclination of 35° and position angle of 125°, which are the values adopted in this study. 

The infrared emission from dust in the diffuse interstellar medium of our Galaxy has 
been extensively studied from IRAS to the Herschel Space Observatory. The quantitative 
Galactic model of iDraine fc Lil (120071 ) consists of a mixture of amorphous silicate grains and 
carbonaceous grains, each with a wide size distribution: sizes range from molecules containing 
tens of atoms to large grains greater than 1 yum. Dust grain temperatures a t the surfaces of 
giant molecular clouds in our Galaxy are hotter than in their cold interiors (IGreenberg fc Li 
19961 ). varying from ~ 15 K to ~ 7 K, respectively. Cold dust grains are a crucial component 
in the ISM of galaxies (e.g.. Block et al. 1994). With reg a rd to the dust temperature - 
metallicity relations in external systems, lEngelbracht et al.l (|2008[ ) find an anticorrelation 
between dust temperature and metallicity. Equilibrium dust temperatures are ~ 23 K for 
solar-type metallicities, increasing to 40 K at a low metallicity of 12 + log{0/H) ~ 8 (see 
their Figure 5). First insights from a Herschel study of Messier 33, with its strong metallicity 
gradient, show that dust temperatures are ~ 25 K in t he central portion of the stellar disk, 
declining to ~ 15 K in the outer domains (Figure 2 in lBraine et al.ll2010l ). 



Here we study the power spectrum (PS) of FIR dust emission from the LMC. Power 
spectra of emission maps are an important but poorly understood diagnostic for interstellar 
structures and the motions that cause them. PS of HI emission from sections of the Milky 
Way are approximately power laws throughout the entire range of observed spatia l frequen- 



cies (ICrovisier &: Dickey 



19831 : lGreenlll993l : iDickev et al 



slopes occur f or PS of Milky Way CO emission (IStiitzki et al. 



/ im emission (IGautier et al. 



1992 



( jDeshpande. Dwarakanath. fc Goss 



2001 



Khalil et al.l 120061 ). Similar 



1998D, IRAS and DIRBE 100 



Schlegel. Finkbeiner. fc DavisI 1998) and HI absorpti 



cion 



20001). These power laws resembl e theoretical expecta- 



tions from passive gas motions (j Goldman! 120001 ) and compressions (ILazarian fc Pogosvan 
2000l ) in a turbulent fluid, so it is natural to interpret them in this way (lElmegreen fc Scalo 
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2004 ). The situation is much more complex for entire galaxies, however, as galaxies un- 
dergo rotation and contain supernov ae, spiral waves, and othe r physical processes that are 
not present in localized clouds. Still, IStanimirovic et al.l (Il999[ ) found a power law covering 
2.5 orders of magnitude in the PS of HI emi ssion from the Small Ma gellanic Cloud (SMC). 
Subsequent analysis of SMC dust emission ( IStanimirovic et al.ll2000l ) gave a similar result. 
In all of these cases, the two-dimensional PS is approximately a single power law with a 
slope between —2. 8 and —3, depending on the tracer and the depth of the velocity slice (e.g., 
Dickev et allboOlh . 



The PS of HI emission from the Large Magellanic Cloud (LMC) is significantly different 
from these others. lElmegreen. Kim, fc Staveley-SmithI (120011 ) found a two-component power- 
law with a break at a wavenumber that is consistent with the inverse of the disk thickness. 
They proposed that turbulent processes and other gas motions are approximately 2D on 
scales larger than the disk thickness, and similar processes are 3D on scales smaller than the 
disk thickness. PS for integrated emission from turbulent gas has a slope that differs by 1 in 
these two cases, with the steeper slope correspond ing to 3D motions at large wavenumbers 
(ILazarian fc PogosyanI l200d : iLazarian et al.ll200ll ). In the LMC, the break scale increases 
with radiu s, indicative of a rea sonable level of flare in the HI disk. Subsequent models of 
disk PS by iPadoan et al.l (120011 ) confirmed this interpretation. 



The same two-part PS was also found for discrete clouds in the Milky Way. iMiville-Deschenes et al. 



(120031 ) found it in HI emission from the Ursa Major Galactic cirrus and modeled it with a 
fractal Brownian motion cloud. Their model suggested that the break in the PS occurs at a 
wavenumber k ^ 1/(2D ) for c loud depth D, and that the power-law slope changes by 1 at 



the break. Ilngalls et al.l ( 12004 ) observed a two-component PS using Spitzer IR observations 
of a molecular cloud in the GUM nebula. They found that the PS slope changed from —2.6 
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at wavenumbers less than 2 x 10 arcsec to —3.5 at wavenumbers greater than 4 x 10 
arcsec"^. The inferred cloud depth was 0.3 pc. 

A break in the PS of HI emi ssion has not generally been observed in galaxies other 
than the LMC. iDutta et al.l (120081 ) studied HI emission from the face-on galaxy NGC 628 
and derived a PS slope of —1.7 without a break. Their PS covered o nly a factor of 10 in 
scale, which may be too small to resolve out two separate components. iDutta et al.l (l2009al ) 
saw a break in NGC 1058 with a steepening of the slope from —1 to —2.5 at an inverse 
wavenumber of 1.5 kpc, which they converted to a disk thickness of 490 pc. Their PS also 
spanned a factor of 10 in scale, but the galaxy may be too distant (10 Mpc) to resolve 
the disk thickness with HI (their angular resolution was 13.6", corresponding to 660 pc). 
PS of Ha and HI emission from several dwarf galaxies s howed power laws and no breaks 
Jwillett. Elmegreen. fc Hunteil I2OO5I : ISegum. et allbood loutta et al.l [2009b[ ) . In the SMC 



-4- 



there is no break either. For the SMC, this lac k of a break could be because the line-of-sight 
depth is comparable to the transverse length (IWesterlundl 119971 ) . This could be the case for 
the dwarf galaxies also (IRoychowdhury et al.ll2010[ ). 



PS miss the break if the high frequency part is dominated by the point spread func- 
tion of bright unresolved emission s ources. These sour ces introduce a Gaussian-shaped dip 
at high frequency in th e PS fe.g.. iDickev et al.l 120011 ). Models of this dip are shown in 
Elme green, et al. ( 2003 ) and Block et al. ( 20091 ). To see a double power law, the part at 
high frequency has to extend for a long wavenumber range from the inverse thickness of the 
disk to the inverse spatial scale of the telescope resolution. If the thickness is smaller than 
the scale of the resolution, then the high frequency power law from 3D turbulence cannot 
be observed. This problem is present for conventional ~ 10" HI resolution in galaxies at 
distances greater than 2 Mpc. Our Spitzer infrared observations of the LMC have a resolu- 
tion scale at least 10 times smaller than the disk thickness, so we should be able to see any 
high-frequency power law clearly. 

Optical observations of galaxy PS have sufficient resolution to see the thickness, but 
unless there are large numbers of stellar groupings with sizes spanning a wide range around 
the disk thickness, the clarity of the PS break will be poor. A tenta tive PS break was 
obse rved optically in NGC 5055 (lElmegreen. Elmegreen &: Leitnerl |2003| ) . but not in NGC 
628 (lElmegreen et al.ll2006l ). A break in the slope of the autocorrelation function for young 
stars i n several nearby galaxies was inferred to mark a change in geometric properties of the 
gas bv lOdekonI bo08( ]. 



Spitzer Space Telescope observations have such high angular resolution that 3.5-4 or- 
ders of magnitude in spatial scale can be covered for the LMC. This is much larger than 
is possible for any other galaxy in any other gas tracer at the present time. The smallest 
resolvable scale is only a few parsecs, which is much smaller than the disk thickness, so there 
is no contamination in the PS from the point spread function. Thus the FIR observations 
of the LMC are id eal for an unpreceden ted view of correlated structures in a galaxy. In a 
companion paper (IBournaud et al.ll2010l ). we describe self-gravitating hydrodynamic simu- 
lations of galaxies like the LMC, with and without energy input from star formation. We 
obtain a two component PS for projected density structure in the simulations too, and we 
examine the velocities and energy sources that lead to this structure. 



In what follows, the observations are presented in Section 12. H the PS are shown in 
Section \'2.2\ while Section [3] contains a discussion of the implications of our results. The 
conclusions are in Section |H 
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Observations 



2.1. Data 



O ur data for the LMC come from SAGE ("Surveying tlie Agents of a Galaxy's Evolu- 
tion"; iMeixner et al.ll2006l ). The images at 24, 70 and 160 fim were obtained with MIPS on 
the Spitzer Space Telescope; they cover ~ 8 x 8 degrees. The total integration time is 217 
hours, which improves the point source sensitivity over previous surveys (e.g., IRAS) by 3 
orders of magnitude. In Spitzer, all bands sample the point-spread function at better than 
the Nyquist frequency (0.4 X/D at 24 and 160 /im, 0.3 X/D at 70 /im). Thus, the instrument 
allows the telescop e-limited resolution of 6", 18", and 40" at 24, 70, and 160 /im, respectively 
(IRieke et al.l 120041 ). The first of these, 6", corresponds to a spatial scale of 1.45 pc in the 
LMC. The theoretical Nyquist limit is 2 pixels, but in reality, the true Nyquist limit is 2 
resolution elements. 

At 24 fim and 70 /im, the images sizes we use are 8192 x 8192 pixels at a scale of 4"98 
and 4"8 per pixel, respectively (the SAGE 24 /xm image was binned 2 x 2 to limit computer 
memory requirements). At 160 /im, the image size is 2048 x2048 pixels at a scale of 15"6/px. 
As far as deprojecting images of the LMC are concerned, there is no unique "center" - the 
dynamic al center of the H I is offset by almost a full degree from the photometric center of 
the bar (IWesterlundl 119971 ) . We choose to conduct our deprojections about the dynamical 
centre of the LMC, which has right ascension a = 5''27.6™' and declination 6 = — 69°52' 



(J2000.0), following section 7 in Ivan der Marel et all (l2002f l. 



2.2. Power Spectra 

Figure [1] shows the MIPS image at 160 fim and the 2D PS at the bottom right. The PS 
is the sum of the squares of the real and imaginary parts of the 2D Fourier transform (made 
with f ourn from Numerical Recipes), averaged over each 2D wavenumber k = {k1 + k'^^y/'^. 
We divided the PS into 15 wavenumber intervals and plotted the arithmetic mean for each 
interval as a dot in Figure [H The PS has 2 distinct parts with different power law slopes. 
The steep slope at high k is —3.08 ± 0.13 and the shallow slope at low k is —2.15 ± 0.48. 
The break is at approximately 1/k = 200 parsecs. Double power laws also occur in the PS 
of the other MIPS channels (Figs. [2] and [3D. The break in the middle of the PS on a log 
scale allows both power laws to be seen clearly; there are 1.5 — 2 orders of magnitude of 
wavenumber on each side for the shorter wavelength bands. 

Emission at 160 /xm arises from cold dust in dark clouds. At 70 /xm, the emission 
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comes from both warm and cold dust, and at 24 /im the emission is from warm dust and 
PAH's. The figures indicate that all of this dust has correlated structure spanning the entire 
resolvable disk of the LMC The PS break is at 100 — 200 pc for all 3 MIPS bands, suggesting 
that the disk thickness is about this value for each dust component. The least-square slopes 
are similar but with a slight progression. At high spatial frequency, they are —3.08 ± 0.13, 
—2.97 ±0.52 and —2.55 ±0.34 for 160 /im, 70 fim and 24 /xm, respectively, and at low spatial 
frequency, they are —2.15 ± 0.48, —1.83 ± 0.36 and —0.78 ± 0.19 for these three passbands. 
The steepening of the slope with FIR passband indicates that there is relatively more spectral 
power on large scales for longer wavelengths. This suggests that cool dust is more diffuse 
and less structured on small scales than warm dust. This conclusion is true regardless of the 
total fluxes and column densities of cool and warm dust, as these quantities contribute only 
to the k = component of the PS. 

For comparison, the 2D PS of HI emission from the LMC has a slope of ~ —3.7 at high 



spati al frequency and a slope of ~ —2.7 at low spatial frequency ( lElmegreen. Kim, fc Staveley-Smith 



200ll ). The break is at a scale of 1/A; ~ 100 pc, increasing slightly with galactocentric radius. 
These slopes are both higher than the corresponding slopes for the dust emission. The HI 
PS is most like the 160 /im PS, although even steeper, suggesting that HI is even more 
dominated by large scale structure than cool dust. 

In order to test the robustness of our analysis, we also deprojected the MIPS images 
about the bar center as opposed to the dynamical center of the LMC; the resulting PS at 24 
Hm remained unchanged, while the slopes of the two power laws at 70 fim and 160 fim were 
only marginally affected. We conclude that choice of center for deprojection purposes is not 
responsible for the presence of two distinct power laws. 



Discussion 



The figures indicate that the relative thickness of the LMC dust disk compared to its 
diameter is only a few percent. This value comes from the ratio of the smallest wavenumber 
in the PS (the inverse of the largest scale) to the wavenumber at the break in the PS (the 
inverse of the disk thickness). The smallness of this ratio indicates that the LMC geometry 
is effectively 2D for some processes. 

Turbulence in a strictly 2D medium (infinitely thin) differs from turbulence in a 3D 
medium. For incompressible turbulence, an infinitely thin layer can have an inverse cascade of 
energy, i.e., energy moves from smaller to lar ger scales, and a direct cascad e of mean-square d 
vorticity, i.e., from larger to smaller scales ( iKraichnanl 119671 ; iLeithI Il968l : iBatchelorl Il969l ) . 
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For 3D turbulence, the energy generally cascades to smaller scales. Thus 2D turbulence on 
large scales can be fed by small scale motions, while 3D turbulence on large scales can only 
be fed by larger scale motions (in the absence of a magnetic field or other physical processes 
that force a connection between large and small scales). 

It would be interesting if an inverse cascade of turbulence occurred in a very thin galaxy. 
Then localized energy from gravitational instabilities on the scale of the Jeans length, which 
is the disk scale height, and from supernovae, superbubbles, and so on, could power large 
scale motions and density correlations. A galactic disk is not strictly 2D even if it is thin, 
but the gas motions can still be highly anisotropic, making the disk effectively 2D. For 
example, density-wave and bar-driven streaming motions are often 5 to 10 times faster 
than the perpendicular motions that produce the disk thickness. Large scale motions may 



be ge nerated by disk self-gravity, t idal forces from companion galaxies (iMastropietro et al. 



2009 ). intergalactic ram pressure ( Tonnesen fc Bryan 2009 : Dutta et al. boiol ). and other 
forcings. They contribute to large-scale structure and to the low-A; part of the density PS. 
More localized energy input from stars and OB associations should d rive 3D motions and 
struc tures that contribute to the high-/c part of the density PS (see also lHodge fc Deshpande 
20061 ). Turbulence driven on large scales would have a shallow PS because of the 2D nature 
of the resulting flows, while turbulence driven on small scales would have a steep PS because 
of the 3D nature of its associated flows. 



Bournaud et al.l (120101 ) ran simulations of galaxies like the LMC and found a two- 



component PS similar to what we observe here. This result occurred whether or not there 
was input from supernova energy, suggesting that gravitational energy alone can drive galac- 
tic turbulence over a wide range of scales. The study also suggested that there may be some 
contribution to the large-scale motions from an inverse cascade of energy in 2D. For example, 
they found a PS for the mean-squa red vorticity, as in nu merical experiments of strictly 
2D turbulence in protoplanetary disks (jPeterson et al.l 120071 ). Thus the LMC PS could be 
the result of gravity-driven turbulence and other large-scale motions contributing to gas 
structure in an effectively 2D medium wider than the disk thickness, supernovae, cloud-scale 
gravity and other small-scale motions contributing to turbulence in a 3D medium within the 
disk thickness, and an inverse cascade of energy at the thickness scale up to larger scales, 
contributing to turbulence throughout the disk. 



4. Conclusions 



The LMC is the only disk galaxy where we can currently study interstellar motions on 
a range of scales covering a factor of 10^ or more. An essential component of this study is 
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the Spitzer Space Telescope, which provides the highest possible angular resolution for large- 
scale maps of interstellar emission. Using the SAGE images at 24 /im, 70 /im, and 160 /xm, 
we derived 2D power spectra of dust emission over the whole LMC disk. The results showed 
one power law covering a factor of 100 in scales at low spatial frequency and another, steeper, 
power law covering another factor of 100 in scales at high frequency. The low frequency power 
law had been observed before with HI emission, but the high frequency power law was more 



limited before, spanning only a factor of ~ 10 in range (lElmegreen. Kim, fc Staveley-Smith 



200ll ). The current observations confirm the two-component nature of the interstellar mo- 
tions, suggesting that gas motions and turbulence are two-dimensional on scales larger than 
the disk thickness, and three-dimensional on scales smaller than the disk thickness. The 
line-of-sight thickness is measured at the break point, and is between 100 pc and 200 pc, 
depending on the wavelength of observation. 
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Fig. 1. — Left: The LMC at 160 fim (North is up, East is left). The brightest part of this 
image has an intensity of ~ 1500 MJy sr~^. The bright cloud in the lower left is the giant 
molecular and atomic cloud south of 30 Doradus. Upper Right: Two-dimensional Fourier 
transform of the deprojected image. The rectangle betrays low level striping in the image 
on the left. Lower Right: 2D PS with power law fits cx A;" and cx at high and low k, 
respectively; the slopes are indicated at the bottom left. Dots with error bars (±lo") are 
averages over a range in log/c. The units of the PS are (MJysr~^)^ for Spitzer calibrated 
images. The smallest scale in the PS is 7.6 pc, which is 2 pixels (31.2") and the theoretical 
Nyquist limit. The PS break at ~ 200 pc (top scale) is interpreted as the line-of-sight depth 
of the dust disk in the LMC. 



Fourier Transform 




LMC - 70 microns 
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10^ 1000 100 10 
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Fig. 2. — Same as Figure 1, but at 70 /xm. The brightest region has an intensity of 4000 
MJy sr~^^. The smallest scale in the PS is 2.3 pc, which is 2 pixels (9.6") and the theoretical 
Nyquist limit. The power spectrum spans 3.5 orders of magnitude. 
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lO 
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Relative Spatial Frequency i^/k^^^) 

Fig. 3. — Image and PS of the LMC at 24 /iin. The brightest region has an intensity of 2800 
MJy sr"^. The smallest scale in the PS is 2.4 pc, which is 2 pixels (10.0") and the theoretical 
Nyquist limit. 



